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We use moving mesh hydrodynamical simulations to make maps of Sunyaev-Zeldovich eflect. 
We present these maps for several cosmological models and explore their lowest moments. We 
find that the first moment, the mean Compton y parameter, is typically between 1 — 2 x fO"'' for 
cluster abundance normalized models, the lower value corresponding to the high density models, 
and scales as ffg"*. Rms fluctuations at 10' scale have an amplitude AT/T 1 — 3 x lO"'^ in 
the Rayleigh- Jeans regime. The amplitude of the power spectrum strongly depends on the power 
spectrum normalization and scales as o-g~^. On smaller scales (/ > 1000) the spectrum is dominated 
by halos below IO^^Mq and so is sensitive to the thermal history of galactic halo gas. On larger 
scales (/ < 1000) the power spectrum is less sensitive to nongravitational energy injection, but 
becomes very noisy, with large field to field variations in the power spectrum caused by rare bright 
sources in the maps, which dominate in the spectrum over the large scale structure correlations. 
Cross-correlation power spectrum with weak lensing or projected galaxy distribution is significant 
and the cross- correlation coefficient is around 0.5 over a wide range of scales. Comparison with 
Press-Schechter predictions gives very good agreement for all the statistics, indicating that there is 
no significant contribution to SZ from non-virialized structures. The one point distribution function 
shows clear deviations from gaussianity and can be well-approximated as log-normal on small scales. 



I. INTRODUCTION 



Cosmic microwave background (CMB) photons propagating through the universe are scattered by hot electrons 
along their path. The effect of this scattering on the photon distribution was first described by Sunyaev and Zeldovich 
and is called the Sunyaev-Zeldovich (SZ) effect. The effect conserves the number of photons, but changes their 
energy distribution. Low frequency photons in the Rayleigh- Jeans regime gain energy on average and are moved into 
the high frequency part of the Planck distribution, with a zero crossing at 217GHz in the nonrelativistic case. The 
amplitude of the effect is proportional to the product of electron temperature and density, or pressure, and is thus 
dominated by hot dense structures such as clusters. 

The SZ effect is nowadays routinely detected by a number of instruments such as BIMA, Diabolo, OVRO, 
PRONAOS, Ryle, SEST, SuZie and Viper (see refs. [|,| for review). All of these have so far concentrated their 
efforts on known clusters, but future experiments will also observe blank fields of the sky with an effort to isolate this 
effect. Among these are the proposed one square degree survey Q and a number of CMB experiments, for which SZ 
may be a significant source of fluctuations that needs to be separated from the primary signal. The Planck satellite for 
example, with its combination of frequency coverage, all-sky map and angular resolution, should be able to identify 
several thousand SZ sources 

The SZ effect has a number of potentially observable consequences in addition to the one descibed above on the 
known clusters. First, the mean SZ distortion can be observed through the deviation of the photon spectrum from 
the Planck distribution. No such distortion has been detected by the COBE/FIRAS data yielding the upper limit on 
Comptonization parameter y (defined below) of y < 1.5 x 10~^ 1^. This is a constraint that has to be satisfied by 
any viable cosmological model. As we will show below cluster abundance normalized models explored here typically 
satisfy this constraint. 

Second, most of the CMB experiments measure fluctuations, for which the most relevant quantity is their power 
spectrum. Again no clear detection of SZ has been reported so far, although this is not unexpected given that the 
main source of fluctuations on degree scales probed so far are primary CMB fluctuations generated at the time of 
recombination. The next generation of CMB experiments is beginning to detect fluctuations at smaller angular scales. 
With higher sensitivities and broader frequency coverage SZ may become detectable even in random patches of the 
sky surveyed by these experiments. There are two related questions that one would like to address. First is the overall 
level of SZ fluctuations and whether SZ may be a significant source of CMB foreground contamination. Second is 
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the power of SZ to distinguish between cosmological models using power spectrum or any other statistics, such as 
counting of sources or nongaussian signatures. 

There are two main theoretical approaches that can be used when calculating the SZ effect. One is analytical 
and is based on the Press-Schechter approximation. In this approach all the mass in the universe is distributed into 
virialized halos. Gas in these halos is assumed to follow a prescribed distribution in density and temperature (usually 
isothermal with a /3 density profile). One then integrates over the mass function and over the redshifts to calculate 
the first moment of the distribution (the y parameter). A similar calculation also yields the uncorrelated contribution 
to the second moment to which one must also add the correlations between the halos . This approach involves 
a number of approximations that must be verified with numerical simulations. For example, some fraction of the gas 
may not be in virialized halos, but may be residing in more diffuse filaments or at the outskirts of the halos. There may 
be significant temperature and density fluctuations of the gas inside the halos. Furthermore, N-body simulations that 
have examined the accuracy of the Press-Schechter mass function have shown it to be a reasonable approximation for 
massive halos, but to overpredict the number of halos at lower masses. The resulting spectrum may be very sensitive 
to the details of the mass function distribution. 

The alternative approach adopted in this paper is to use numerical simulations to calculate the SZ effect. Our 
approach is based on moving mesh hydrodynamical simulation code JTo[ and we use ray tracing through the simulation 
box to generate maps of SZ. This is the approach that was adopted previously by [ |ll|Jl^ . Most of these early studies 
were of low resolution, but nevertheless showed that SZ could be an important effect on the CMB. Recently 
generated SZ maps using SPH simulations of higher resolution, although they did not calculate the CMB power 
spectra. An alternative method based on 3 — power spectra was developed in Q and was recently adopted to the 
MMH code as used in this paper by [Esl. We compare our results to the previous work where possible. 



II. METHOD 

We use moving mesh hydro (MMH) code developed by one of us The code is a hybrid between Eulerian 

and Lagrangian grid based hydrodynamic methods. By deforming the grid in the dense regions along potential flow 
lines it provides a ten fold increase in resolution compared to fixed grid Eulerian codes, while maintaining regular 
grid conditions everywhere pO| ] . The code has been succesfully parallelized on shared memory systems and has a low 
computational cost per grid cell. A 256'^ run takes around 2-4 days of wall clock time on 32 processors. 

We ran about a dozen simulations varying cosmological parameters, box sizes and particle/mesh dimensions to 
check for various numerical effects. Our main models are tCDM, KGDM and OCDM. Their parameters are 

• tCDM: flat model with = 1, f^A = 0, cts = 0.56, VLbh = 0.034 and = 0.25. 

• ACDM: flat model with f7„ = 0.37, f^A ^ 0.63, as = 0.8, ^th ^ 0.034 and n^h = 0.25. 

• OCDM: open model with a„ = 0.37, Ha = 0.0, ag = 0.8, ^bh = 0.034 and a„,h = 0.25. 

We ran several simulations varying box sizes from 50ft.~"'^Mpc to 200/i^^Mpc to verify the effects of mass and force 
resolution on small scales and lack of power on large scales. We also varied mass and scale resolution, using mesh 
sizes between 128^ and 256"^ and typically placing 1-8 dark matter particles per mesh cell. Our largest simulations 
are ACDM and OCDM models with 100/i"^Mpc box, 256^ mesh and 256^ dark matter particles. 

During the simulation we store 2-d projections through the 3-d box at every conformal time step that corresponds 
to a light crossing time through the box. The projections are made alternatively in x, y and z direction to minimize 
the repetition of the same structures in projection. We store projections of SZ, kinetic SZ, gas and dark matter 
density. For SZ we store the projection 

Ay = JL^n.T^a^x 1 

meC^ 

where ut is the Thomson cross-section, c the speed of light, TOg the electron mass. Ax the comoving width of the box, 
a the expansion factor, rig electron number density and Tg electron temperature, both of which are obtained from 
the output of the MMH hydro code (note that we ignore relativistic corrections, which are only relevant for rare hot 
clusters). Our 2-d maps are 1024^ for 256'^ and 512^ for 128'^. We have verified that this preserves all the information 
content by comparing the results to the higher resolution projection. Number of stacked projections depends on the 
box size and ranges between 25 (for tCDM with 200/i"^Mpc box) to 150 (for ACDM with 50/i"^Mpc box size). 

After the simulation is completed we use the 2-d projections to make maps of SZ. We stack together the maps of 
SZ separated by the width of simulation box, randomly choosing the center of each box (note that use of periodic 
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boundary conditions guarantees there is no edge for any of the maps). We then project these maps onto a map 
constant in angular size. In principle the angular scale of the projection should be determined by the angular scale 
of the simulation box at the initial z, since any larger scale produces repetition of the same structures in the map. 
However, most of the structures in the maps are coming from low z and these cover very little volume of the simulation 
box; hence they do not repeat itself even if we increase the angular scale beyond the size of the box at the highest 
redshift. Typically we project up to z ^ 4 — 30 using oversampling between 1-4 without any noticeable artifacts in 
the maps. We checked for possible artifacts in the power spectra by comparing weak lensing maps produced in the 
same way to the analytic predictions [ p^lJl^ . The agreement was very good in all models. For SZ it should be even 
better, since SZ power spectrum is much more dominated by isolated sources and, as will be shown below, large scale 
correlations do not dominate even on large scales. 

A given simulation can be used not only for the simulated parameters, but also for models with lower fig. To do 
this one can simply relabel an earlier time output as being one at a later time by changing the redshift of projection 

..,ncw /-^old I 1 \ .....new /^old -i 

where superscript old stands for original simulation and new for the new one. Once the redshifts are relabeled one can 
project over the 2-d projections using the nearest projection to the required z. For a flat model with no cosmological 
constant this scaling preserves all cosmological parameters (except erg). This scaling is particularly useful to obtain 
the dependence on tJg in a flat model independent of other parameters. We verified its accuracy by running two 
simulations: one with erg = 1 but using the rescaling of equation ^ to obtain erg = 0.56 model, and another with 
(Tg = 0.56 at the final output, and found very good agreement between the two. The simulation with erg = 1 was then 
used to predict cg dependence of the y parameter and SZ power spectrum. 

It is also easy to verify the effect of the baryon density in the limit Jlfc <^ f2„i. In this limit the baryons are just 
a tracer of dark matter potential and do not dynamically couple to it. This fails in the cores of the halos, which 
are however already not correctly modelled in these adiabatic simulations without cooling and non-thermal energy 
injection. In this approximation SZ scales linearly with baryon density times Hubble constant ^hh. We only show 
results that are normalized to BBN nuclesynthesis constraint Of,/i^ = 0.018, but results for other values of rif,/i can be 
easily obtained using this scaling. Incidentally, one could also use transformation of scale to map a given model into 
another family of models, this time varying the shape and the amplitude of the input power spectrum. By combining 
this transformation with the time transformation above one can find a family of models which are normalized to erg 
today and have a different shape of the power spectrum. The tra de-off in this case is a loss in force and mass reso lution. 



Some of the maps produced by the simulations can be found at http: / /fcynman.princcton.cdu/^uros/sz.html 



III. PRESS-SCHECHTER PREDICTIONS 



The Compton parameter y — J dy is given as a projection over the electron pressure along the line of sight (equation 
|l|). This can be reexpressed as a line of sight integral over the density weighted temperature of electrons, 

y = — / « (1 + ^gWx, (3) 



where Sg is the gas overdensity and rfg mean electron density today. 

In Press-Schechter (PS) picture all matter in the universe is divided into halos of a given mass. The mass distribution 
is specified by the halo density mass function dn{M) / dM . This can be written as 

^dM^A/(M)dM, (4) 

where p is the mean matter density of the universe. The function f{M) denotes the fraction of mass in halos of mass 
M. It can be expressed in units in which it has a universal form independent of power spectrum or redshift if we 
express it as a function of peak height v = [Sc{z) / <t{M)]'^ , where Sc is the value of a spherical overdensity at which 
it collapses at z (Sc = 1.68 for Einstein-de Sitter model) and (t(M) is the rms fluctuation in spheres that contain on 
average mass M at initial time, extrapolated using linear theory to z . For scale free spectra 

^, , dn d\nM , , 

The actual form given by Press & Schechter is vf{v) = (i'/27r)-'^/^e~''/^, although modified versions of this form 
that fit better N-body simulations have been proposed ]l9| . Direct hydrodynamic simulations have shown good 
agreement with our form ^ |^ . 
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The virial temperature of the halo in the spherical collapse model is only a function of virial mass and is given 
by ksTe = qM^/^. The conversion factor from mass to temperature q is approximately unity if mass is expressed in 
units of \Q^^h~^MQ and T in keV. Density weighted temperature is given by 

{{l + 5)T)=q j f{v)dvM^/\ (6) 

This gives ((1 + 5)T) = O.SkeV today in a KCDM model with Sl^ = 0.37 and (7% — 0.8, with other models giving 
comparable numbers. Further assuming that gas traces dark matter, ie 5g — 6, we can calculate the integrated 
Compton y parameter from equation |[ This gives y = 2 x 10~^ for this model. Flat models give a factor of 2 lower 
value because of a more rapid evolution of clusters with z. This reduces y relative to a low density model assuming 
the cluster abundance today is the same. 

Using equation]^ we may also explore the dependence of y on as- For this we need to relate a{M) to M. Linear 
power spectrum at the cluster scale can be approximated as a power law P{k) oc cr|fc". Since rms variance tT^(M) 

4 

scales as k^P{k) one finds k oc j/(ct|)~^/^"+^\ Mass goes as M oc pR^ oc k~^ for which we find M^/'^ oc fc~^ oc erg"''" . 
On cluster scales we have —2 <n < —1, hence {{1 + 5)T) oc (Tg~^. The y dependence on a% depends on the projection 
of ((1 + 5)T) along the line of sight (equation but gives qualitatively similar result. 

One can also compute power spectrum using the PS model. To do this one can first calculate the 3-d power 
spectrum of density weighted temperature P(ij^s)T{k) and then project it along the line of sight using the Limber's 
equation. 



Ci = 327r3 



P(i+S)T{k ^l/r,a)-^, (7) 



where r is the comoving angular distance to x, given by Xi i?sinh(x/i?) or _Rsin(x/i?) in a flat, open or closed universe 
with curvature R, respectively. The above expression is valid in Raylcigh- Jeans limit where AT/T = j{x)y — —2y. 
This can easily be rescaled to another frequency using the spectral function j{x) — x{e^ + l)(e^ — 1)^^ ^ 4 with 
X = hv/kBT (which in the limit a; ^ gives j = —2). 

In the PS formalism there are two contributions to the pressure power spectrum [Q-^ 



P{k)^P'^{k) + P'"'{k) (8) 

The first term pP{k) arises from the correlations within the single halo. This term contribution to the 3-d power 
spectrum is given by 



Pil+S)Tik) = ^ / f{iy)d:.-M^/'\y{k)\\ (9) 



{2ny J p 

where y{k) is the Fourier transform of the halo profile normalized to unity on large scales (fc — > 0), assumed here for 
simplicity to be independent of Af . One power of M above is given by mass pair weighting and M'^/'^ is given by the 
square of the temperature of the halo. This term is heavily weighted toward rare massive clusters, resulting in a very 
strong Poisson term compared to the halo-halo correlation term discussed below. As shown in this term dominates 
over the halo-halo correlation even on large scales, where it behaves as a white noise. In the large scale limit [k — > 0, 

y{k) = 1) the scaling with erg is given by pP^g^^ oc (7^^/^'+"^ oc al-^\ This will again be modified somewhat by 
the integration over the redshift, but does provide qualitative understaing for the strong cg dependence seen in the 
simulations described below. 

Second term is the contribution from halos correlated with one another. On large scales these cluster according to 
the linear power spectrum Pii,i(fc), except that they can be biased relative to the dark matter. The halo bias can be 
either larger than unity for halos more massive than the nonlinear mass or less than unity for those below that. An 



expression that fits N-body simulations reasonably well was given in |§,|2l[| (see also [g^ 22 1 for a modification relevant 
for lower mass halos) 

&M = 1 + ^- (10) 

Oc 

The halo-halo contribution to the density weighted T power spectrum is given by 



-f(l+(5)T(^) — -Plin(fc) 



q / f{v)di,M^'^{v)b{v)y{k) 



^puk)[{b)(,+s)T{[i+m]\ (11) 
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where (&) (1+5)7- is the pressure weighted bias, defined as 



(b) 



{1+S)T 



J f{v)dvM^''^{v)h{v)y{k) 



(12) 



For KCDM model we find in the large scale limit {y = 1) (6}(i+5)t — 3, hence SZ halos are significantly biased relative 
to the dark matter. This is because the T oc M^/'^ weighting weights preferentially towards the large halo masses 
which are biased. The bias decreases towards smaller scales where y{k) suppression for finite k is more important for 
larger more massive halos and the dominant contribution shifts towards smaller, unbiased or antibiased halos. Despite 
this large bias the halo correlation term is small relative to the halo term even on large scales and Press-Schechter 
model predicts that SZ power spectrum does not trace large scale structure on large scales. 

If SZ does not trace LSS, the large scale bias b cannot be measured from its power spectrum. Cross-correlating SZ 
with weak lensing or galaxy map is a more promising way to obtain this large scale bias. This is because the Poisson 
term is given by M^^^ weighting. 



Ml 



(1+5)T,5 



(fc) 



(27r)3 



/(j,)d,.^M2/3|y(fc)|^ 
P 



(13) 



and is so less dominated by rare objects than SZ power spectrum with M^/'^ weighting. The correlated contribution 
to the cross-correlation between density weighted T and dark matter density is given by 



f{v)dyM'^'^b{M) 



f{v)dyb{M) 



^PUk){b)ii+s)T{{l + S)T), 



(14) 



where the second integral above is the mass weighted bias, which is by definition unity. Numerical comparison between 
the two terms indeed shows that the are comparable on large scales. To obtain the SZ-weak lensing or galaxy cross- 
correlation power spectrum we again use this 3-d power spectrum in Limber's equation ^ with the appropriate window 
function for the dark matter or galaxy projection. 



IV. SIMULATION RESULTS 



A. Mean Compton parameter 

We first present the results on the mean y parameter. The value for this parameter is 1.0 x 10^^ for tCDM in 
lOO/i-^Mpc simulations with 128^ mesh, 2.4 x 10"^ for OCDM in lOO/j-^Mpc with 256^ mesh and 2.0 x 10"^ for 
KCDM with lOO/i-^Mpc with 256^ mesh or SO/i^^Mpc with 128^ mesh. For the latter y decreases by roughly 20% if 
the resolution is decreased to 128"^ mesh, suggesting that small halos not resolved in larger box simulations contribute 
a significant fraction to its value. This is also confirmed by Press-Schechter type calculations [Q. Non-thermal energy 
injection and cooling may provide additional corrections to the results of the adiabatic simulations. For example, 
p3[ include feedback in their simulations and find significantly higher density averaged temperature (IkeV) than our 
simulations (0.3keV), although other simulations that also include heating do not (G. Bryan, private communication). 
The mean y parameter is therefore sensitive to the thermal history of the gas. 

Our results are still one order of magnitude below the current experimental limits from COBE/FIRAS y < 1.5 x 10"^ 
, so it is unlikely that this constraint will play a major role in distinguishing between the models that differ in the 
history of energy injection in the universe. On the other hand mean y parameter does put a constraint on the ag. 
This scaling of y with can be obtained from the same simulation using the method described in previous section. 
We find y oc o-\~^ , which is in a good agreement with the Press-Schechter predictions. To violate the FIRAS limits 
one needs to increase cg by a factor of 2 over the cluster abundance value. COBE normalized standard CDM with 
(Tg — 1.3 would be problematic. Our results are a factor of 2 lower from those in for the two low density models 
(once we account for the difference in Qbh), while we are in a good agreement for tCDM model. This disagreement 
does not seem to be explained by limited mass resolution on small scales, since our simulations have higher mass and 
force resolution, yet predict lower values for y for the low density models. We agree well with results in ||l5[| . 
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B. SZ power spectrum 



Next we explore the power spectra of SZ in the Rayleigh- Jeans regime. We first investigate the effects of resolution. 
These are shown in figure |l| for different parameters of the simulation for KCDM model. Shown are 256^ 100/i~^Mpc 
simulation, 128^ 100/i~^Mpc simulation, 128^ 200/i~^Mpc simulation and 128^ 50/i~^Mpc simulation. They aU agree 
on large scales, where the shape is close to the white noise model, characteristic of a power spectrum dominated by 
rare sources. On small scales 200/i~^ simulation begins to loose power for I > 1000 compared to other simulations. 
Comparison between 128^ and 256^ 100/i^^Mpc shows that the two are in good agreement on scales up to Z ~ 1500, 
beyond which 128'^ begins to loose power. Similar conclusion is obtained by comparing 50/i~^Mpc and 100ft."^Mpc 
both with 128^. The figure shows that a 50/i^^Mpc box is sufficient for the power spectrum on large scales. This 
indicates that the large scale power spectrum is not dominated by very massive clusters above IO^^Mq which do 
not form in such small boxes, but rather by less massive IG^^^^^Mq clusters, which happen to be nearby in the 
projection along the line of sight. This conclusion is also confirmed by the Press-Schechter calculation, where the 
brightest sources correspond to this mass range. On small scales the 50/i~^Mpc box has a power spectrum comparable 
to 256^ 100/i~^Mpc simulation with the same mass and force resolution. Very little small scale power therefore arises 
from mode-mode coupling with scales larger than 50ft.~^Mpc. Overall 100/i~^Mpc 256^ simulation resolves the power 
spectrum between 100 < I < 5000, while 128^ simulations of the same size are sufficient between 100 < I < 2000. 
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FIG. 1. Power spectrum comparison between different resolution and box size simulations for KCDM model. Also shown 
is the white noise spectrum scaling as 1^ . On large scales the SZ spectrum is white noise in all cases. 



On large scales the power spectrum approaches white noise in slope. This is an indication that the power spectrum 
is dominated by the rare bright sources in the map and not by the correlations between them, which would give a much 
shallower slope. The same behaviour is seen also in other cosmological models, shown in figure ||. While OCDM gives 
similar predictions to KCDM model, tCDM is significantly lower, just like in the case of mean y parameter. Low 
density KCDM predicts 3 x 10~^ fluctuations around / ^ 100, which rises to 3xl0~^ at / ^ 5000. Comparison with 
primary CMB in figure ^ shows that the SZ is unlikely to contaminate CMB power spectrum for MAP, while Planck 
and smaller scale experiments should be able to measure the SZ power spectrum. Our results are in good agreement 
with lis], although we note that our OCDM spectrum is somewhat higher than theirs, in better agreement with the 
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PS calculations. 
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FIG. 2. Power spectrum comparison between cluster normalized tCDM, XCDM and OCDM. The latter two have 
significantly more power. Also shown is the primary CMB (long dashed) which dominates on large scales. 

Figure ^ shows field to field variations in the power spectrum for 2° x 2° maps. It shows significant fluctuations on 
large scales. This is another indication that the power spectrum on large scales is dominated by rare bright sources. 
On smaller scales the power spectrum is more robust. The main contribution on those scales is from smaller more 
abundant halos, so there is less sampling variance. Note that as the map size increases the chance of finding a bright 
source in it increases as well, so the median power spectrum on large scales grows with the map size. This is equivalent 
to the power spectrum after the brightest sources have been removed. The power spectrum in a typical few degree 
map should correspond to the power spectrum with a few thousand brightest sources removed. As shown in this 
can reduce the power spectrum on large scales significantly and should be kept in mind when comparing our results 
to other predictions. For example, power spectra in are based on all the power without the removal of bright 
sources, which should in general give a higher amplitude on large scales. 
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FIG. 3. Power spectrum variations for several 2° x 2° fields in XCDM model. There are a factor of 2 variations from field 
to field even on large scales, caused by rare bright objects in the map. 

Figure ^ shows the change in the power spectrum as a function of erg for tCDM model. The amplitude of the power 
spectrum is very sensitive to as. Doubling this parameter changes the power spectrum by 2 orders of magnitude. 
Fitting to a power law we find C; oc Cg. This steep dependence on a% can be understood with Press-Schechter 
formalism developed in §3 and is caused by a rapid increase in number of bright sources as a function of trg: and is 
in contradiction to the estimates of jl2j . One can also see a crude estimate of the scaling relation by noting that the 
Ci oc< T > 6"^ and the temperature is a function of the non-linear length scale tnl, i-e. T oc (rNL^^o)^ fUl- For 
CDM like spectra near the non-linear mass scale we have tnl Cg"^, giving Ci oc cTg"^' similar to what we find 
in simulations and to the PS arguments. This indicates that care must be exercised when extracting cosmological 
parameters from SZ maps, since similar differences arise between low and high density cosmological models (figure |^). 
A small increase in a% of the order of 20% changes the spectrum as much as does changing the density from Vl„i ~ 1 
to Q,m = 0.3. Unless we are confident that we know the local value of cg to better than this accuracy we cannot use 
SZ to infer the density of the universe. Similar argument also explains the redshift dependence of SZ power spectrum, 
which is very strong. Most of the contribution to the power spectrum comes from z < 1. This is also true in OCDM 
and ACDM, although the z dependence is less steep there. 
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FIG. 4. Dependence of SZ power spectrum on as in tCDM model. 



C. Non-Gaussian signatures 

Since we have 2-d maps of SZ we may use these to study non-Gaussian signatures in SZ. There are two reasons why 
to investigate non-Gaussian signatures of SZ. First is that such information can be used to determine the cosmological 
model. Example of this is the SZ luminosity function, where one identifies SZ sources and plots the number density 
as a function of the flux in SZ . The abundance of such sources in the limit where they are unresolved has been 
made by the Press-Schechter formalism or peak-patch formalism |^^. It has been shown that the number 

density of sources is strongly sensitive to the density parameter, similar to the mean y and power spectrum statistics 
discussed above. However, the sources have some internal structure, they may be clustered, may contain substructure 
and may not be spherical, all of which complicates such analytic approaches and they need to be verified using the 
simulations. This is of particular interest to the proposed surveys of small fields centered on a random portion of the 
sky which will not in general detect massive clusters, but rather a collection of smaller sources. Simulations such 
as these presented here are necessary for investigation of the non-Gaussian effects in such random portions of the sky. 

Second reason to study non-Gaussian signatures is that they may provide additional leverage in separating the SZ 
from the primary CMB, at least under the assumption that primary CMB is Gaussian. Even when one does not have 
sufficient frequency coverage to distinguish the two components on the basis of their different frequency dependence 
one can use non-Gaussian signature of SZ to estimate its contribution to the power spectrum. The only assumption 
in this procedure is that at a given smoothing scale there is a strong correlation between the non-Gaussian signature 
and its second moment, so that the latter can be estimated from the former. 

It was shown above that the 2-point statistics is dominated by rare bright sources and so is very noisy on large 
scales. This will be even worse if one considers higher order statistics such as skewness and kurtosis. For this reason 
we concentrate here on the one point distribution function (pdf) of SZ smoothed at a given angular scale. These 
are shown in figure || for a 2° x 2° field in ACDM model with several smoothing radii, all smoothed with a top-hat 
window. As shown in figure || the pdf is approximately lognormal on small scales, in the sense that there is an excess 
of large y decrements caused by bright rare sources. The pdf can be transformed into an approximate Gaussian 
if plotted against log(j/). For larger smoothing angles the pdf becomes narrower and approaches a Gaussian, so it 
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becomes more difficult to distinguish it from primary CMB on this basis. 




FIG. 5. The one point distribution function of logy for several smoothing radii with a top-hat window. 



D. SZ-weak lensing and galaxy cross-correlation 



SZ map can be cross-correlated with other maps, such as weak lensing, galaxy or X-ray maps. The first two are 
sensitive to projected density and may not be so strongly dominated by rare bright objects, while the latter traces 
projected p^T^I"^ and is even more strongly dominated by rare objects than SZ. Examples of weak lensing maps are 
those reconstructed from shape distortions of z ~ 1 galaxies or from CMB distortions a,t z ^ 1100. Example of 
projected galaxy maps are those from APM, SDSS or from numerous degree field surveys. They can be parametrized 
by the mean galaxy redshift. In the case of SDSS one can use photometric information to weight the galaxies according 
to their distance to optimize the signal to noise of the cross-correlation. Even more promising are smaller and deeper 
surveys, such as those used for weak lensing studies with several hundred thousand galaxies over a degree area. 

A useful quantity to compare the cross-correlation between different maps is to compute the cross-correlation 
coefficient 

Figure ^ shows the cross-correlation coefficient for weak lensing map reconstructed from background galaxies at 
z — 0.5,1,1100. The cross-correlation coefficients for all 3 cases is between 0.4 and 0.6 and only weakly depends 
on scale. Similar results are also obtained if one cross-correlates SZ with a galaxy catalog with a mean redshift of 
z ^ 0.2 — 0.5. This result shows that both SZ and weak lensing maps are dominated by relatively nearby objects, 
which is why the cross-correlation coefficient does not significantly decrease as the redshift of background galaxies is 
decreased. This is good news for shallow surveys such as SDSS, which should be able to detect the cross-correlation 
signal when compared to MAP or Planck maps. 
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FIG. 6. Cross-correlation coefficient between SZ and weak lensing convergence assuming background source is at 
z = 0.5, 1, 1100. Similar results are obtained for cross-correlation with galaxies whose mean distance is at half the distance to 
the lens source. 

The amplitude of the cross-correlations can also be compared to the analytic predictions developed in §3 and we 
find good agreement between the two. The predictions have also been made in the case of weak lensing of CMB by 
1^,^. Our results are several times lower on scales between 100 < I < 3000. A more detailed comparison shows that 
the assumed bias in the two models is comparable (6 ~ 3 — 4) and the difference is mostly caused by lower density 
weighted temperature in our simulations (0.3keV), consistent with PS, compared to their assumption (IkeV). Density 
weighted temperature is sensitive to the thermal injection from stars and supernovae and it could be significantly 
higher than in our simulations, which neglect this effect. However, such heating of the gas would shift the weight in SZ 
to cooler halos (below IkeV), which are not biased or are even antibiased. So non-thermal increase in density weighted 
temperature also requires a decrease in bias and the two effects nearly cancel out, leading to a smaller senitivity of the 
cross-correlation power spectrum on the density weighted temperature than one would naively assume. Unfortunately 
this also implies that the expected signal to noise of SZ and CMB cross-correlation will be significantly lower than 
predicted [^6|j2j] and becomes only marginal using our results. 




Using hydrodynamical simulations we have produced maps of SZ and analyzed some of their low order statistics. 
The mean Compton parameter was found to be in the range y ~ 1 — 2.5 x 10~^, an order of magnitude below current 
FIRAS limits. Only if one increases as by a factor of 2 over the cluster abundance constraint does one violate the 
FIRAS limit, something which is clearly excluded based on cluster abundance data. FIRAS limits will therefore not 
play a major role in constraining the cosmological models. 

The power spectra were found to be noisy and white noise like on large scales, indicating that they are dominated 




V. DISCUSSION AND CONCLUSIONS 
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by uncorrelated bright sources. The amphtude at I ~ 1000 is 1 — 3 x 10~^ and does not represent a major source of 
foreground to primary CMB on large scales {I < 1000). On smaller scales SZ power spectrum should be detectectable 
and should dominate over CMB for I > 2000. On these scales the predictions become sensitive to the thermal energy 
injection into the gas and simulations that ignore such effects become unreliable. Conversely, power spectrum of SZ 
on small scales should give us important information on the thermal history of the gas in small halos. Note that 
the dominant contribution to the power spectrum comes from relatively nearby structures with z < 1, so SZ will not 
provide information on gas history at high z. 

On scales where SZ is not negligible compared to primary CMB one can use the non-Gaussian signatures of SZ to 
estimate its contribution to the power spectrum. We have shown that on arcminute scales the one point distribution 
function is well approximated as a log-normal, differing significantly from a normal distribution. This signature can 
be used to estimate the contrubution of SZ to the power spectrum even in the absence of multifrequency information. 

Further insight into the thermal history of the gas can be obtained by cross-correlating SZ with other maps that 
trace large scale structure, such as weak lensing or projected galaxy map. This can also provide information on 
correlations between groups and clusters on large scales and their bias relative to dark matter. The cross-correlation 
coefficient is quite high, of order 0.5, across a wide range of scales even when the redshift distribution of galaxies or 
lensing mass peaks well below z ^ 1. This strong correlation provides further incentive to planned SZ surveys on 
random patches of the sky, since even if sensitivity is not sufficient to detect SZ directly, one may be able to detect it 
through cross-correlation with weak lensing or galaxy maps, thus providing information on the thermal state of the 
gas in the universe. 

We thank Eiichiro Komatsu, Alexandre Refregier and David Spergel for useful correspondence. U.S. ackowledges the 
support of NASA grant NAG5-8084. Computing support from the National Center for Supercomputing Applications 
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